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DIGITAL LINE-STORE STANDARDS CONVERSION: PRELIIVllfSlARY INTERPOLATION STUDY 



Summary 

The experimental Digital Line-Store Standards Converter^ built in Research 
Department served to demonstrate the good picture quality obtainable with digital 
techniques, even when using a comparatively simple linear interpolation aperture utilising 
information from only two lines of the input picture. The work described in this report 
was a study, mainly theoretical, into the requirements for a more sophisticated interpola- 
tion system involving the use of signal contributions from more than two lines. 
Computer optimisation techniques were used in determining the most suitable aperture 
shape using contributions from different numbers of successive field lines, and a block 
schematic diagram was derived for a proposed experimental interpolator to be used with 
the experimental digital line-store converter. A series of subjective tests is in hand to 
establish the basic interpolation parameters of a converter for service use. The results of 
these tests will be given in a further report. 



1. Introduction 

The use of digital techniques in standards conversion 
offers the potential advantages of high reliability, small 
size, reduced cost, and, provided that the arithmetical 
operations involved in signal processing are correctly carried 
out, optimum picture quality from the output signal, how- 
ever complex the signal processing. This last advantage is 
of particular importance when considering the requirements 
for interpolation, since, as a rule, the more lines which are 
available simultaneously, the more accurate will be the 
estimate that can be made of the picture information that 
would have been produced had the signal been originally 
scanned using the output standard. 

Information from a block of lines can be made 
available simultaneously by using the signals appearing at 
the inputs and outputs of a number of one-line delays in 
series; analogue signal delays of this duration produce small 
but perceptible picture defects. However, shift registers 
provide a means of delaying digital signals with no degrada- 
tion and it is possible to achieve fairly sophisticated inter- 
polation without involving undue cost and complexity. 

The equipment proposed in this report is not intended 
as a design prototype for future digital line-store converters, 
but rather as a research tool to enable subjective tests to be 
carried out with the purpose of determining the practical 
limitations which can be imposed on the complexity of the 



interpolation whilst achieving a high standard of output 
picture quality. 



2. Theoretical considerations 

2.1. A definition of interpolation aperture 

Interpolation is necessary in standards conversion 
because the scanning of the original scene samples the 
picture vertically at a rate (in terms of spatial frequency) 
which depends on the number of lines per field. If a tele- 
vision picture is to be obtained as if generated on another 
standard, the original scene must first be recovered from the 
input samples to effectively give a continuous function 
which can then be re-sampled, and hence rescanned. 

According to Nyquist's sampling theorem, if a signal 
is sampled at a rate greater than twice the highest frequency 
it contains, it can, in theory, be recovered perfectly from 
the samples using a suitable filter. However, it can, in 
practice, be recovered with various degrees of accuracy by 
replacing each sample, represented as an impulse, by an 
aperture function whose origin coincides with the impulse 
position and whose magnitude is scaled in proportion to the 
strength of the impulse. The value of the recovered signal 
at any point is then the sum of the various aperture 
functions corresponding to the train of samples. A repre- 
sentation of this process is shown in Fig. 1 where the 
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original signal recovered 
by interpolation from samples 



samples 



E 2 =A 2 Si etc: 

where A^is ttie value of the aperture 
J] [r. function corresponding to the distance R,Q. 




~*\ distance between p— 
successive scanning lines 

Fig. 1 - Interpolation as a summation of displaced aperture functions. The recovered signal 

at Q is the sum of E^, E^, E^ and E^ 



original signal has been recovered by interpolation from 
the samples. It will be observed that the aperture function 
used here as an example resembles (sin x)/x over a limited 
range. The choice of a suitable practical aperture will be 
discussed in Section 2.5. 

The aperture function need not be symmetrical about 
its origin; but if it is, an alternative approach is to consider 
it centred on the point at which the interpolated value is 
required. The magnitudes of the neighbouring impulses 
are then scaled by the relative value of the aperture at the 
various impulse positions and the scaled values summed to 
produce the required interpolated value as shown in Fig. 2. 



The process of displacing and weighting a function is a 
convolution operation and thus interpolation can also be 
described as multiplication in the frequency domain where 
the spectrum of the sampled signal is simply multiplied by 
the spectrum of the interpolation aperture. The action of 
sampling a signal causes its spectrum to appear as sidebands 
centred on multiples of the sampling frequency as in Fig. 3. 
The spectra centred on the sampling frequency and its 
harmonics can, however, be almost completely rejected if 
the spectrum of the interpolation aperture approximates the 
pass-band of a low-pass filter that would transmit only the 
zero-order spectrum representing the original signal. 



£^ = Ai S^ 

£2 = A2 S2 
etc 




Fig. 2 - Interpolation as a summation of weighted samples. The disturbance at Q is equal to the sum of E , E , E and E 
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Fig. 3 - Spectrum of signal sampled at a frequency f 
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Fig. 4 - Alias components generated when original signal contains spectral components 

at frequencies greater than f^/2 
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F/gr. 5 - 777e spectral behaviour of sampled, interpolated 
and resampled signals 

(a) Original signal (h) Sampled signal 

(c) interlace-sampled signal (d) and (e) Result of interpolating 

ib) and {c] with an 'ideal' aperture 

(/) and {g) Resampled and interlace-resampled interpolated signals 

(h) Addition of resampled and interlace-resampled signals 

(/) Resampled signal subjected to restricted bandwidth) 



2.2. Space-time duality 

The signal representing tine brightness down a vertical 
line in the scene has spatial frequency components and is 
spatially sampled. Thus distance can be measured in terms 
of picture height (p.h.) or line height |l.h.) and frequency 
in cycles/p.h. or cycles/l.h. However, there is a complete 
analogy with the time domain in that, for a particular 
element position along a scanning line, the samples occur at 
regular intervals of time (line period) and the interpolation 
aperture becomes the temporal impulse response of a low- 
pass filter. This is the approach used in the design of 
analogue line-store converters now in service. 

2.3. The idea! aperture 

From Fig. 3 the ideal interpolation aperture would 
seem to be that which has the rectangular spectrum of an 
ideal low-pass filter cutting at half the sampling frequency, 
but other factors may influence the choice. 

(1) The original scene may contain spectral components 
at frequencies greater than half the sampling fre- 
quency, and these components will be accentuated if 
excessive vertical aperture correction has been em- 
ployed. The process illustrated in Fig. 3 will then 
produce overlapping of the repeated spectra, an effect 
known as 'aliasing' since the spectral components 
overlaid are not distinguishable from the signal com- 
ponents. In such cases it may be desirable to reduce 
the cut-off frequency of the aperture spectrum to 
eliminate the alias components as in Fig. 4. 

(2) The vertical bandwidth will be too high if inter- 
polation is followed by down-conversion. In this 
case it may be desirable that the cut-off should be 
reduced below half the input sampling frequency. 

(3) The effect of interlace. 

This third factor is by far the most important and is dis- 
cussed further. 

The following argument applies to down-conversion 
when both input and output standards are interlaced. 

Suppose that the spatial spectrum of the brightness 
signal down a vertical line in the original scene is as shown 
in (a) of Fig. 5. If the scanning action of the source at the 
input standard takes place with a spatial frequency /| c/p.h. 
assuming that the source aperture is infinitely small the 
spectrum of the sampled signal is shown in (b); as can be 
seen, higher order spectra are produced. The spectrum of 



the second, interlaced, field has inverted components (as 
explained in Reference 4) for odd orders as shown in (c). 
If spectra {b) and (c) are combined to form a complete 
picture at the input standard the even orders reinforce each 
other and the odd orders partially cancel each other out 
to an extent which depends on the temporal response of the 
display phosphor. Hence, integrated over a complete 
picture the effective value of /, is doubled. However, 
because the eye is able to resolve temporally the individual 
fields, the partial cancellation of the odd orders gives rise 
to the familiar phenomenon of 'strobing' at low input 
frequencies (e.g. a plain raster) and 'twittering' at high 
frequencies (e.g. edges). 

Suppose, however, that spectra (b) and (c) are inter- 
polated and resampled in an interlaced manner before being 
combined. If the interpolation aperture has a rectangular 
spectrum cutting at/^, where/^</j/2 the spectra of the inter- 
polated fields are shown in (d) and (e). 

If the interpolated signals are now resampled with a 
spatial frequency /q, the spectra of the two fields are shown 
in (f) and (§■). When spectra (/I and t^) are combined to 
form the complete output picture as in (h) steady and 
alternating components are produced as outlined above. 

Finally, a limit must be imposed on the range of dis- 
played frequencies because of visual acuity or the existence 
of line broadening caused by the finite display spot size. 
A reasonable assumption is that these effects restrict the 
spectral content of the displayed picture to an upper 
frequency of /^ as shown in {/). 

The final result can be illustrated by examining the 
output for a sinusoidal input whose frequency varies from 
zero to /|. This can be done by tracing through the 
products derived from a particular spectral component in 
Fig. 5(a). For an input frequency / in the range zero to 
/j. indicated by F^ in Fig. 5(a) an alternating component at 
/q— / is produced in addition to the steady component at/. 
If/ lies between f^ and f—f^ no output is possible but if/ 
lies beyond f—f^ as indicated by F^ in Fig. 6(a) a steady 
and an alternating spurious component are again produced. 
As the dotted lower curves in Fig. 5(/) are parts of the curve 
of Fig. 5(a) translated by an amount/ and the solid lower 
curves are parts of the same curve translated by an amount 
(/'j-/q) the alternating component is at a frequency /-/ 
and the steady component is at a frequency /— (/j— /q)- 
This behaviour is summarised in Fig. 6. The alternating 
component at f^-f in region (1) would be present in a 
signal generated at the output standard and is therefore 
inherent. If, in region 3, the steady spurious component is 
the more objectionable and its visibility decreases with 
increasing frequency then there might be an advantage in 
lowering/,,, thereby narrowing region 3 and delaying the 
onset of the steady component. This would be paid for by 
a loss in definition. On balance there does not seem to be 
any strong reason for not putting /j, equal to//2. It must 
be emphasised, however, that although it is possibly the 
best solution it is not ideal because of the spurious products 
in region 3 and therefore the result of the conversion process 
field by field can never be the same between interlaced 
systems as that obtained by generation at the output 
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standard at the source. The aperture profile having this 
flat spectrum cutting off at/j/2 is of (smx)lx form, having 
zeros spaced at the distance between successive lines of an 
input field. 

2.4. Limitations on the practical aperture 

If the best aperture is of (sin;ic)/x form it is clear that 
a practical aperture of limited width and having an imperfect 
spectrum approaches the 'ideal' as its width increases. 
There are, however, certain conditions of varying impor- 
tance which should preferably be satisfied by any aperture 
of limited width. These conditions are discussed below. 

(i) The aperture spectrum must have zeros at integer 
multiples of the input sampling frequency /. If this 
is so the higher orders of a sampled constant signal 
(corresponding to a picture having uniform vertical 
brightness) are completely rejected and no spurious 
components are produced. Failure to meet this con- 
dition would, after resampling, give rise to horizontal 
striations on a uniform background. This condition 
is also equivalent to the requirement that values of the 
aperture profile spaced apart a distance of one field 
line should sum to a constant regardless of their 
position, an extension of Equation 2 in Reference 4. 

(ii) The value of the constant in condition (i) should be 
unity if the interpolation process has unity gain. 

(iii) The spectrum of the aperture should have zeros at 
integral multiples of the output sampling frequency so 
that no spurious product is formed at zero frequency 
except when the input frequency is/. In general, a 
spurious component at zero frequency can only occur 
by the translation in the resampling process after 
interpolation of a frequency at an integral multiple of 
f^. Thus if the aperture spectrum has zeros at these 
frequencies, no spurious dc component can result. 



(iv) The value of the normalised spectrum at half the input 
sampling frequency should approximate to 0'5, corres- 
ponding to the ideal case. This point is discussed 
further in Appendix A. 

(v) The aperture profile should have no discontinuity at 
the ends. See, for example. Fig. 2. Failure to meet 
this condition with apertures having a width not equal 
to an integral number of lines implies a discontinuity 
in output signal as the position of the output line 
varies, assuming a uniform d.c. input signal. With 
apertures having a width equal to an integral number 
of lines no discontinuity of output signal occurs, but 
an extra sample is needed to produce the output 
when an output line coincides with an input line. 

(vi) The gradient of the profile should be zero at the ends 
(except for any discontinuity) if it is continuous else- 
where (Fig. 2). Failure to meet this condition with 
apertures of both integral and non-integral line widths 
implies a discontinuity of gradient at points on the 
profile spaced at integral line widths from the ends. 



2.5. The derivation of a practical aperture 

There are several possible approaches to the problem 
of deciding on the profile of a practical aperture. The 
approach described here is a limited Fourier series method, 
an extension of a method already described elsewhere.'' 



where w is the aperture width in p.h. and d is the distance 
coordinate in p.h. measured relative to the centre of the 
aperture. 

Fig. 7 is drawn for an aperture width equal to an 
integral number of lines but the method of building up the 
spectrum applies equally well to non-integral line widths. 
When the width is an integral number of lines, so that 



W=(« + 1)//j 

where n is the order of interpolation. 



(3) 



then the spectrum has the value ^„ + .j/2 at the frequency 
/j which can thus be specified as zero thereby satisfying 
condition (i) of Section 2.4. 

Moreover, to simplify the mathematics, the values of 
A^ for r>«+1 can also be set to zero thereby ensuring that 
the spectrum passes through zero at the points rfJ{n+'\) 
for jr|>«+1. Thus the aperture profile becomes a limited 
Fourier series given by 

S{e)=A^ +A^ cose +A.:^ cos 20 +A^ cosnO 



n 

= ^^^cosre,|6'|<7r 
r= 



and the aperture spectrum F{p) is given by 



(4) 



F{p) =A^ sine vw +—i 



sine iv — ^/w)w + sine (j' + Viv)w H- 



The starting point is a theorem^ which states that an 
aperture of width w units has a spectrum the value of which 
can only be independently specified at discrete frequencies 
which are multiples of ^/w cycles/unit. The spectrum in 
between these frequencies is derived by the superposition 
of appropriately scaled (sin x)/x functions centred on each 
frequency, and having zeros ^/w cycles/unit apart. This 
follows because the aperture can be expressed as the pro- 
duct of a repetitive waveform having a period of w units 
and a rectangular pulse of width w units and the spectrum 
of the aperture is the convolution of the spectra of these 
two functions. This process is illustrated in Fig. 7 for an 
aperture width w equal to 6//| p.h. where 1//| is the input 
sampling spatial period. The spectra of the repetitive 
waveform and the pulse are shown in {a} and {b) respec- 
tively and the spectrum of the aperture in {d). As can be 
seen the spectrum parameters^Q,^^ etc. are the amplitudes 
of the Fourier components of the repetitive waveform, at 
the frequencies 0, ^^w etc. Thus the aperture profile can 
be expressed as the Fourier series 



r= 



sine {v - "/w)w + sine [v + 



sine (pw - r) + sine (vw +r)\ 



"/w)w 



(5) 



where sine vw = sin(7rpw)/(7rwj 
rectangular pulse of width w. 



is the spectrum of a 



S(d) =A^+A^ cos6i +A^cos2d ..., \6\<n 



(1) 



in Equation (1) 6 is an angular coordinate which varies 
from — 7r to +tx across the aperture and is given by 



d = 2-n d/w 



(2) 



2.6. Optimisation of the aperture 

The profile defined by Equation (4) satisfies condition 
(i) of Section 2.4 because ^„ + ^ is zero and condition (vi) 
because it is a finite series. However it does not necessarily 
satisfy the remaining four conditions as can be seen by 
referring to Fig. 8 which shows an aperture of width 4 lines 
and the problem of optimisation is to decide on the impor- 
tance of the various conditions consistent with obtaining 
a spectrum as close as possible to the 'ideal' rectangular 
form of cutting at/j/2. If there are insufficient degrees of 
freedom because n is too small it will be impossible to 
impose all the conditions simultaneously. 
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F/jj. 7 - Derivation of aperture spectrum of widtfi 6/f-^ p.fi. 

(a) Spectrum of repetitive waveform (b) Spectrum of single square pulse of width 6/f-. 

(c) Some components of aperture spectrum (dj Aperture spectrum obtained by summation of all components 



Assuming «>2 so that all the conditions can be 
accommodated there are four linear equations relating 

Aq Ajj and the spectrum is optimised as a function of 

the remaining independent parameters. 

The optimisation can be performed by minimising an 
error function which represents the deviation of the actual 
spectrum from the ideal spectrum. An error function which 
has been used in theoretical work was formed by summing 
values of the spectrum expressed in decibels at regular 
frequency intervals within the passband (i.e. up to/j/2) and 
adding to this a similar contribution from the stopband but 
expressed as a deviation from a desired value of —40 dB. 



Alternatively the parameters A^ A^^ can be 

given arbitrary values and the spectrum and profile in- 
spected to see how closely the four conditions are obeyed. 
If necessary the parameters can then be adjusted. This 
second approach was found to be more successful particu- 
larly as it was found in practical work that some of the 
conditions could be relaxed. 



2.7. Spatial quantisation of the aperture 

The ratio of 625 to 405 line frequencies is such that 
81 output lines occur in the same time interval as 125 input 
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Fig. 8 - Typical aperture spectrum and profile obtained by optimisation process 
{a) Spectrum (b) Profile 



distance l.h. 



lines. Tills relationship can also be assumed to hold for the 
spatial positions of the lines. This means that only 81 
discrete values of the interpolation aperture profile per line 
of each input field need be known. If, however, the 
aperture is several lines wide the number of values that 
needs to be stored may become prohibitively large and, 
moreover, if the values are to be set independently as part 
of a subjective experimental arrangement the complexity of 
the apparatus reaches impractical proportions. 

The required number of stored values can be reduced 
if certain restrictions are imposed on the interpolation 
aperture. Firstly the aperture can be assumed symmetrical 
so that the number of required values is halved. Secondly, 
the apertures can be spatially quantised into less than 81 
discrete values although impairments in interpolation will 
then result in addition to those caused by the finite width 
of the aperture. The form of these impairments can 
easily be analysed in the frequency domain. 

The quantising can be expressed as a sampling action 
followed by a boxcar operation (Fig. 9). Thus, if the 
aperture is quantised into k quanta per input line height it 
is equivalent to sampling with a frequency/ = kf-^ c/p.h. 
and forming each sample into a rectangular pulse of width 
Ikf^ p.h. In spectral terms if the original aperture spectrum 
is F{v) the quantisation produces a spectrum 



F^{v) = %mc vikf- ^ F(v + tkf{, 



(6) 



where sincM = sin (■nu)l{-nu). 
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Fig. 9 - Derivation of spatially quantised aperture for k = 3 
(a) Original aperture (b) Sampled aperture (c) Quantiseci aperture 
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Fig. 10 - Derivation of ttie spectrum of a spatially quantised aperture fork = 3 
{a) Original spectrum (b) Spectrum of sampled aperture (c) Spectrum of quantised aperture 



The summation term results from the samphng, and the 
sine term from the boxcar operation. The derivation of a 
typical form of F Av) is shown in Fig. 10. As can be seen 
the effect of quantisation is to introduce additional spectral 
lobes at the frequencies ±kf-^, ±2kf^, etc. 

Now the first scanning process (at the input frequency 
/j) produces components centred on 0, ±f^, ±2/j etc (Fig. 5) 
and thus significant energy exists near the frequencies ±fc/|, 
2fc/j . , .''. after interpolation. The effective resampling 
process described in Section 2.3 (at the output frequency 
/q] translates all components through integral multiples of 
/g and thus the energy in the additional lobes may even- 
tually produce substantial components near the origin. It 
is shown in Appendix B that the result of the two 
sampling processes, with or without intervening interpo- 
lation, produces additional spurious components at spatial 
frequencies separated from the input frequency by integral 
multiples of /|/125{= /q/SI). If the corriponent at 

/;•(= / + ''/i/125) is termed the rth resultant order the 
amplitude of this component, e(f^), is given by 



e(/;) = 



Fq(f+[35/-- 



81i] /:) 



(7) 



This equation also holds for unquantised apertures 
when Fq becomes F but because of the absence of the 
additional lobes, the level of the spurious products is much 
lower. It is also shown in Appendix B that Equation (7) 
reduces to the unquantised case when k = 81. This must be 
so if interpolation requires values of the aperture only at 
the centre of each quantum. This argument can be 
extended to cover the case when k is any multiple of 81. 



The energy pattern of resultant orders repeats every 
81 so that the 81 orders nearest the origin contain all the 
information about the levels of the spurious products. The 
effect of quantisation can now be assessed by comparing 
the level of these 81 products with and without quan- 
tisation. 

T he physical significance of the resultant orders is best 
appreciated by considering an input signal in the form of 
sloping bars of sinusoidal cross-section for example the 
corner bars of a Test Card. It is shown in Appendix B that 
the spurious components also appear as sloping bars as 
shown in Fig. 11. 

The interaction of the spurious bars with the input 
bars gives rise to the familiar phenomenon known as 
'knotting'. 

For example, when linear interpolation is used the 
level of the relevant resultant order causing the strongest 
spurious component is about -20 dO and this is visible on 
the corner bars of a Test Card as can be seen in Fig. 12. 

2.8. Theoretical results obtained 

Aperture spectra have been optimised according to 
the criterion mentioned in Section 2.6 for various integral 
values of K up to 7, that is for aperture profiles up to 8 lines 
wide. In general, increasing the value' of n gives a sharper 
cut, as might be expected. For example, with « = 7 the 
spectrum falls from 90% to 10% of the value at the origin in 
a frequency range of 0- 1 1/, compared with 0-32/| for the 
present analogue converters.* 

* This information was communicated to the authors by Mr. R.E. 
Davies. 
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Fig. 11 - Positions of maximum brightness on the display 

for an input of sloping bars showing the —44th order 

spurious component. The output scanning lines are 

not shown 

To obtain theoretically satisfactory rejection of un- 
wanted signals it was found necessary to relax condition (iv) 
of Section 2.4, allowing the spectrum to take on values 
much lower than 0-5 at the frequency /j/2. In the absence 
of condition (iv) however, the optimisation resulted in 
spectra having a very good stopband response but a poor 
passband response. 

The second approach of setting arbitrary parameter 
values was then tried. The imposition of condition (v) was 
found to be most restricting and when this was removed 
spectra were produced which were theoretically and experi- 
mentally acceptable and which generally satisfied condition 
(iii) quite closely. Finally condition (ii) was relaxed very 
slightly to allow for the fact that the aperture values were 
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Fig. 12 - Photograph illustrating interpolation errors in 
the corner bars of Test Card C 



quantised in amplitude and therefore did not always sum 
to unity. Fig. 8 shows a typical aperture spectrum and 
profile obtained. 

As a subsidiary investigation certain profiles were 
quantised into 9 and 3 quanta per line and the resulting 
spurious products calculated and compared with those from 
the unquantised aperture and the quantised ideal aperture. 
A summary of the results is shown in Appendix C where 
components down to the —40 dB level are included. 

it is seen in Appendix C that the behaviour of an 
aperture having a width of 8 lines is virtually identical to 
that of an infinite-width aperture. Moreover, it can be 
shown that for 9 quanta per line the spurious products for 
an infinite-width aperture are at least 29 dB down. As 
products at —20 dB cause an impairment of about grade 3 
(see Section 2.7) it is considered that an aperture profile 
having 9 quanta per line should be acceptable. 

2.9. The effect of contrast law on interpofation 

The above theoretical treatment is based upon one 
important assumption which is not necessarily upheld, 
namely that the chosen interpolation aperture will operate 
upon signals having a linear relationship to brightness. In 
practice, standards converters are located at transmitter 
stations, and so receive an input composite signal which has 
been gamma-corrected. 

Standards conversion of gamma-corrected signals will 
generally result in an incorrect summation of aperture 
samples, these errors being greatest when the difference in 
brightness from line to line is greatest. 

A series of exploratory tests was carried out using the 
digital converter with a linear interpolation aperture and 
suitable switching arrangements were provided whereby 
conversion could take place before or after gamma correc- 
tion, and the pictures produced by the two processing paths 
were viewed alternately. Slides covering a good cross 
section of subjects were selected (including those known to 
show up interpolation errors), and a suitable high quality 
slide scanner was used as a source of signals. However, 
very little or no improvement could be detected when con- 
verting signals before gamma correction. 

It was concluded that for the majority of subjects, 
the line to line differences in picture content were restricted 
by the vertical resolution of the system, and were conse- 
quently never large enough to show up the inaccuracies 
caused by operating on gamma-corrected signals. It was 
provisionally concluded, therefore, that gamma correction 
was unlikely to influence the design of the interpolation 
system, although further confirmatory tests will be carried 
out when more sophisticated interpolation methods are 
available. 



3. Practical proposals 

To assess the practical importance of the theoretical 
conclusions reached above, subjective tests must be per- 



10 



formed on the visible impairments to the television picture 
caused both by simple and also by more refined inter- 
polation methods. For this work a flexible and compre- 
hensive interpolator is being developed which will provide 
immediate comparison of the same programme material 
interpolated in various ways. 

Fig. 13 shows the basic block diagram of the proposed 
interpolator. The tapped delay will consist of seven one- 
line digital delays so that signals from eight adjacent field 
lines of the picture may be simultaneously available for pro- 
cessing. The eight video signals will be multiplied by 
appropriate coefficients in eight multipliers and the resulting 
products added, with due regard to sign, in the adder. 



625/50 digital 



picture signal 



I 1 

tapped delay 



adder 



multipliers 



?*» 



+ 



405/50 
signal out to 



time-redistributor 



V 



interpolation 
change over 



coefficient 
store 1 



coefficient 
store 2 



sync. 1625/50 - 
pulses]405/50- 



coefficient 
selection control 



Fig. 13 - Basic blocl< diagram of proposed interpolator 



By suitable choice of the coefficients loaded into a 
store, the width, shape, amplitude resolution and to some 
extent the spatial quantisation of the aperture may be 
varied. It is hoped that a set of interpolation parameters 
will be found which are both subjectively satisfactory and 
simple to instrument in operational converters. 



4. ConcEusions 

The best aperture spectrum for standards conversion 
between interlaced television systems is probably rectan- 
gular, cutting at half the input standard field spatial fre- 
quency. The corresponding aperture profile is of (sinx)/x 
form having zeros one field line apart. A practical aperture 
can approach this ideal as closely as desired as its width 
increases. If the aperture is expressed as one cycle of a 
limited Fourier series the resulting spectrum is a function of 
a limited number of variables and can be optimised subject 
to certain constraints. These constraints are that the spatial 
spectrum should have zeros at integral multiples of the 
input and output spatial scanning frequencies and fall to 
half amplitude at half the input spatial scanning frequency. 
In addition the aperture profile and its gradient should both 
fall to zero at the ends of the aperture. If the aperture is 
spatially quantised to save storage of information the 
smallest number of quanta per input-field line that gives 
acceptable degradation probably decreases with increasing 
aperture width. For a width of 8 input lines 9 quanta is 
probably acceptable. The minimum width of aperture, 
amplitude resolution and number of spatial quanta per line 
which are acceptable remain to be determined subjectively. 
To this end an experimental interpolator is being construc- 
ted which, it is hoped, will provide a means for exploring 
the subjective effect of varying all aperture parameters. 
The results of these tests, giving recommendations to 
Designs Department, will be the subject of a further report. 



The complexity of the equipment increases rapidly as 
the width of the aperture, the number of bits in the co- 
efficients, and the fineness of the spatial quantising of the 
aperture are increased. For example, an aperture eight 
lines wide, with eight bits per coefficient and nine quantised 
positions, requires an interpolator using nearly 1500 simple 
integrated circuits with correspondingly high cost. For an 
experimental interpolator this degree of complexity has 
therefore been taken as the upper limit. 

Immediate comparison between two chosen sets of 
interpolation parameters will be provided by the inter- 
polation changeover switch which will select coefficients 
from one or other of two coefficient stores in use. A 
number of plug-in store boards will be available by means 
of which comparisons can be quickly made. 
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APPENDIX A 



The Transmission of a Siiarp-Cut Filter at the Cut-Off Frequency and Related Theorems 



Suppose a d.c. signal is sampled so that its spectrum 
consists of lines at integral multiples of the sampling 
frequency, /g. If this signal is passed through a filter having 
a rectangular spectral response cutting at ±f^ where f^<f^ 
then only the zero-order spectral component is transmitted 
and the output is d.c. In the time domain the output of the 
filter can be expressed by convolution as the superposition 
of (sin x)/x pulses of unit height having TJ2 between 
successive zeros and spaced apart a distance T^. If, there- 
fore, Tj = aTj2 = aT^ the output of the filter is given by 



enit) 



n = — oo 



sin 7r(? +naT^)IT^ 
naTJIT 



Esin TTJi 
T!(t + 



= constant for 0<a<2 (8) 



In particular for a = 1 each peak of a pulse coincides with 
the zeros of all the others and therefore e^ at the peaks 
must be unity. Thus 



sin iT{t + nT^)IT^ 
, Ti{t + nT^)/T^ 



(9) 



From energy considerations the value of the constant in 
Equation (8) must be proportional to /^ and therefore 
inversely proportional to a. Thus 



sin IT {t + mT^)/T2 ^ 
irit + naTJ/T, a 



for 0<fl<2 



(10) 



When a = 2 pulses occur at intervals of IT^and conse- 
quently each peak again coincides with zeros. Thus e^ is 

unity at the points 0, 2T^, 4T^ However, e^ is zero 

at the points T^, 37^, because pulse zeros occur at 

these points. Further, e^ must be symmetrical about the 

points 0, 2T^. AT^ and T^, 3T^ because of the 

symmetry of the pulses and therefore the gradient must be 
zero at all these points. Turning to the spectral domain, 
when a = 2,f^ =f^ and consequently the first-order spectral 
lines fall at the cut-off frequency. Together with the 
previous argument it is clear that e^ involves a d.c. and a 
sinusoidal term such that 

e^{t) = 72(1+ cos 2nf^t) (11) 

The corresponding Fourier transform of e^ is 

V2{b{v) + yAb(v-f^)+b(v+f^]]] (12) 

where b(p) is the Dirac-Delta function, showing that the 
spectral lines at ±f^ have half the amplitude of that at the 
origin. Thus, if the spectrum of the filter is given by 



F(v)= 1,|y|</e 
F(v) = Q,\v\>f 



(13) 



the value of F{v) at \v\ = f^ is 0-5. Any real filter which 
approximates to F{v) should, therefore, also approximate 
to the value of 0-5 at the cut-off frequency. 

When 2<a<4 the first-order spectral lines fall inside 
the pass-band and thus the output of the filter is 



e^{t) = constant (1 + 2 cos 2Trt/T^) 



(14) 



For a = 3, T^ = 3T^ and consequently each peak coincides 
with zeros. Thus e^ = 1 at the points 0, 3T^, 6T^ 

Therefore, fora = 3 



E°° s\mT{t + n3TJ/T . 
— ~ — ? = V3 ( 1 -h 2 cos 27r?/3r ) (15) 
ir{t + n3T^)T^ ^ 

n = —00 z' z 



Again from energy considerations the value of the constant 
in Equation (14) must be inversely proportional to a, thus 



smirit+naTJ/T 1 

^ -- a + 2 cos 2tt t/aT^ 



Esm TT[r 
■n{t + 



(t+naT^)IT^ a 



for 2<fl<4 



(16) 



When a = 4 the second-order spectral lines fall at the cut-off 
frequency and therefore a further term must be included in 
the right-hand side, in general, by a similar argument it can 
be shown that for a i^2p,p integral 



sm -nit +naT^)IT^ 1 



Esin -nyi 
"7(7+ 



1 + 2 cos 2-nt/aT + 



(t+naT^)IT^ a 
+ 2 cos 27r 2tlaT^ + -i- 2 cos 2-nmtlaT^] 



1 +2 y^ co%2TmtlaT^) 



n = 1 



(17) 



where m is the integer next below fl/2. 
For a = 2p,p integral 

E°° sm-n{t+naT^)IT^ 
■n{t+naTJ/T^ 
n = —00 z" z 



- ( 1 + 2 V^ cos 2mtlaT^ + cos 2nl2T^ ) (18) 



Equations (17) and (18) are useful in that they enable an 
infinite series to be expressed in terms of a finite series. 
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APPENDIX B 
Spurious Products Resulting from Two Successive Sampling Operations 



B1. Derivation of the products 

Suppose a sinusoidal signal of frequency / is sannpled 
at a frequency /j. The resulting spectrum consists of pairs 
of spectral lines centred on the points «/| (at nf^ ± f) 
where n is the first sampling order. If this sampled signal is 
then resampled at a frequency /^ each component of the 
spectrum is repeated with a period of /^ so that the result- 
ing spectrum consists of pairs of lines centred on the points 
«/j + m/j, where m is the second sampling order. 



Suppose also that 



/o//i =pIq 



(19) 



where p and q are integers without common factor. Also 
let 



/o = P^ /i = lh 



(20) 



so that h = fjq is a conveniently small unit in terms of 
which to measure frequency differences. Then the result- 
ing (output) spectrum consists of pairs of lines centred on 
the points (nq + mp)h = rfjq, say where r is an integer. 

Now it is well known that if p,q have no common 
factor, then integers m^ n^ always exist so that 



moP+n^q = 1 



(21) 



and it is also clear that these integers are not unique, since 
m^ can be replaced by m^ + \q and n^ can be replaced by 
«Q — Xp, where X is any integer. Further, for any value of 
r we can find integers m J, « J such that 



m^p +n^q =r 



In fact, we can take 



OTj =rmQ — sq 



■m^ + sp 



(22) 



(23) 



where X is any integer, and the«jth first sampling order will 
contribute to the rth resultant order if n^ has any value 
satisfying (23). 



If interpolation follows the first sampling process, the 
amplitude of the nth first sampling order is simply 

F{nf;+f)+F{nf;-f) 

where F{v) is the interpolation aperture spectrum at spatial 
frequency p. Hence the rth resultant order consists of lines 
at frequencies r/|/^ ±f or rh ±f, and the total output is 



e{rh±f)= ^ F(n/-±f) 

S = — oo 

where « J is given by (23). 



(24) 



If the aperture is quantised into k quanta per field 
line, F[u) is given by Equation (6) instead of by Equation 
(5) and the expression for the rth resultant order becomes 

oo 

e{rh±J)= ^ |sinc[(«/|±/)/(fc/i)] X 

S = — oo 
oo 

^ F[(«, +kt)f;±fl\ (25) 



t = -oo 



where, as before, 

h =f-Jq «j =nj+ps . 

For 625/405 conversion, we are mainly concerned with the 
particular case 

p = 81; ^=125; m^ = -54; «o = 35 (26) 



B2. Symmetry considerations 

As the modulus of any spectrum is symmetrical about 
the origin and all spectra considered are real 



e{rh ~f)=e{-rh +f) 



(27) 



The left-hand term is the rth order product resulting from 
the spectral line at —/and the right-hand term is the —rth 
order product resulting from the line at +/. Thus all the 
spectral line amplitudes can be obtained by considering the 
products only from the line at +f. Further if 

/^ = rh +f< 

there will be a product of equal amplitude at — /^. Thus 
considering only positive frequencies, products occur at 
the frequencies j/^|. Thus Equation (25) becomes 



<?(l/rl' = Yl 1''"'' ^'"i-^i +f^lW >« 

S = —oo 
oo 

2^ F[(n^ +kt)f.+f]\ 



(28) 



t = -oo 



B3. Special cases 

For A: = 81 



^ F[(«j +A;r)/i+/l = ^ F[(3&- + 8U)/i+/] 

f = _GiO t = — oo 

which is independent of s. Thus 
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1 y] sinc[(«/i+/)/81/i]J X 

Is = — oo ' 



I Y^ Fl{35r + 81t]f;+fli (29) 

(t = -oo ) 

But the first factor of e in Equation (29) is of the form 

y sine (r + «r)/r which 

n = — oo 

by Equation (9) is equal to unity. Thus 

CX3 OO 

e= ^ F[(35r + 8U)/;+/l = \^ F{n/-^+f) 

t = — OO S = — OO 

which is the unquantised case above. 

Similarly it may be shown that the case k = 162 
reduces to the unquantised case. 

The 'Ideal' Aperture 

i^Q(j') = 1 0<|i;|</i/2 

F^(v) = \v\>fJ2 
Forj/|</i/2, 

OO 

^ F„[(«j +A:f)/i+/] =1for«j +A;?=0 



t = -OO 



= otherwise. 



The condition for F^ to be different from zero may be 
written 

t = -(njk) = -(3Brlk) -(81x/;c) (30) 

As t is an integer there are two special cases. 

(1) k = Z, Q, or 27, s has any value, r is a multiple of k. 
Then only the Q,k,2k . . . . ik orders are non-zero and 
are given by 

OO 

e = ^ sine [ (25ik + SU) /j +/] Ikf-^ 



which by Equation (17) can be written 

k *" 

1 +2 

X = 1 



e =- 



81 



1 + 2 ^ cos 27rs(35fA:/j +/)/81/i 



(31) 



(2) /: = 5, 7 or 35, r has any value, s is a multiple of k. 
Then all orders are non-zero and if s = 0, ^, 2A: . . . . 
ik . . . the rth order is given by 

OO 

e= 2^ sinc[(357- + 81/^)/|+/l/fc/i 

i= _oo 

which by Equation (17) can be written 
40 



1 
81 



1 + 2 ^ cos 27rs(35r/i +f)IBW] 



s = 1 



(32) 



B4. The effect of interlace 

When the first and second sampling processes are both 
interlaced the odd orders produced by each process are 
inverted for the interlaced field. Thus for the interlaced 
field an odd second-sampling order generated from an even 
first-sampling order is inverted but an odd second-sampling 
order generated from an odd first-sampling order is not 
inverted. The various alternatives are summarised in Table 



Now considering the relationship 

81m + 125« =r 

the parity is summarised in the truth table. Table II. Com- 
paring Table I and Table 11 it can be seen that the even 
resultant orders are steady whereas the odd resultant orders 
alternate during a complete picture. 



TABLE I 



TABLE II 



m 


n 


Resultant 








S 





E 


A 


E 





A 


E 


E 


S 



m 


n 


r 








E 





E 





E 








E 


E 


E 



where t is the integer next below 81/{2A:). 



Odd 

E Even 

A Alternating 

S Steady 



B5. Genera! observations 

The general expression for e is best solved by com- 
puter. The infinite summations are limited by placing a 
limit on the modulus of the argument of sine and F. In 
practice limits of 10 and 10/| were found satisfactory. 

In general for unquantised apertures the strength of 
the first sampling orders decreases monotonically as n 
increases. Thus the strongest resultant orders are those 
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which contain contributions from first sampling orders near 
the origin, that is, low values of n. For example for « = 1 
the first sampling order is at/ + /| and contributes to the 
resultant order at this frequency, namely the 125th resultant 
order. By similar reasoning, it is possible to predict the 
sequence of resultant orders in descending strength. Be- 
cause the energy pattern of resultant orders repeats every 
81 it is sufficient to quote the order number of lowest 
modulus. The sequence is 

0, ±37, ±7, ±30, ±14, ±23, ±21 

This sequence breaks down for quantised apertures because 
the strength of the first sampling orders does not decrease 
monotonically due to the presence of the additional 
spectral lobes. 



B6. The physical significance of the products 

The analysis performed so far has involved real spectra 
corresponding to horizontal bars of sinusoidal cross-section 
symmetrically disposed about the horizontal axis. Under 
these conditions the presence of the spurious components 
produces an amplitude modulation of the original signal. 
However, the effect of the spurious components is more 
easily appreciated by considering sloping bars of sinusoidal 
cross-section as in Fig. 1 1. 

Without loss of generality the origin can be placed at 
the intersection of a bar peak and an input scanning line 
making the vertical spectrum of the sampled signal down 
the linex = real as before. For an arbitrary value of x say 
Xq the vertical signal suffers a linear shift J'q where 



and flg is the slope of the bars. 



(33) 



If the vertical frequency of the bars is / the linear 
shift imparts a phase shift of -2-nfy^ to the components 
n/i -1- / and +2Trfy^ to the components at n/j - / in the 
sampled signal. If the signal is now resampled, with or 
without interpolation, and a scanning line passes through 
the origin then each component is repeated with a period 
/q as before. Thus the theory already developed holds 
except that all products resulting from the line at -n/ suffer 
a phase shift of —2-nfy^ and vice versa. Moreover, the 
symmetry conditions apply so that considering only positive 
frequencies, products occur at \f + nf^ +mf^\ = [f -i-/-/j/125| 
= If^l with a phase such that if /^>0 the phase is negative 
and vice versa. For/^0 the phase shift of -2nfy„ pro- 
duces a linear shift oi fyjf^. and therefore the slope of the 
spurious product is fyjfpc^ = ajlf^. For/^<0 the phase 
shift of +2'nfy^ produces a linear shift of -fyj\f,.\ =fyjf,. 
as before. Thus for all values of /^ the spurious component 
has a vertical frequency of \f^\ and a slope a given by 



a = a///, = fl„,/(1+7/i/125/) 



(34) 



For example the corner bars of Test Card C have a 
vertical frequency of 80 c/p.h. and a slope of 0-75. As 
/: = 312y2 c/p.h. Equation (34) becomes 



0-75 
1 -^r/32 



(35) 



From Section B5 the strongest spurious component occurs 
for r = 37, -44, -125 etc. The more objectionable effect 
occurs for even values of r which produce steady spurious 
components. Thus taking r = -44, a = -2 which gives a 
pattern having a slope somewhat greater than the orthogonal 
value of -4/3 and having a coarse vertical frequency of 
30 c/p.h. This pattern is visible when linear interpolation 
is used although the level of the -44th resultant order is 
about -20 dB (calculated from the spectrum of the linear 
interpolation aperture). 
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Width 



Infinite 



Siines'^' 



2 lines 
with 405 
zero 



2 lines 
without 
405 
zero 



APPENDIX C 
Strength of Resultant Orders for an Input Frequency of 0-3 f ' ^ * 



Quanta 
per line 



dB relative 
to input 



Oto-10 
-10 to -20 
-20 to -30 
-30 to -40 



Oto-10 
-10 to -20 
-20 to -30 
-30 to -40 



Oto-10 
-10 to -20 
-20 to -30 
-30 to -40 



Oto-10 
-10 to -20 
-20 to -30 
-30 to -40 



81 



resultant order number 



-37, -7 

+37, +7, +30, -30, -14 

+ 14, +23, -23, -21, +21 


+37 



resultant order number 



-9, +9 

+ 18,-18, +27 

-27, +36, -36 



-9, +9 

+ 18,-18, +27 

-27, +36, -36 



(3) 



+37 

-37, -9, -7, +9, +7 

+18, -18, +30, +27, -27 


+37 

-9, +9, +18, -18, +28 
-35, +27, -27, +36, -36 



resultant order number 





+30, 

-30, -21, +21, +9 

-9, +39, -39, -12, +12, -3, +3, +27, -27 

-24, +24, +6, -6, -36, +36, +15,-15, +33, -33 



+30 

-30,-21, +21, +9,-9 

+39, -39, -12, +12, -3, +3, +27, -27 

-24, +24, +6, -6, -36, +36, +15, -15, +33, -33 



+30, -30, +37, -21 
+21, +9, -9, +39, -39 
-12, +12, +7, +18, -18 



+37 

+30, -30, +7, -21, +21 

-14, +9, -9, -23, +39, -39, +16, -12 

+12, +28, +18, -18, -35, -33, +33, -2 



(3) 



NOTES 

1. This frequency gives most spurious components for the 3 quanta case 



2. 

3. 



Each order strength is 0-5 dB greater than the infinite width case 
+37th order also significant at a frequency of 0-4/| 
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